First-principles study of native point defects in topological insulator Bi 2 Se 3 
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The p-type Bi2Se3 is much desirable as a promising thermoelectric material and topological in- 
sulator, while the naturally grown Bi2Se3 is always n-type doped by native point defects. Here we 
use first-principles calculations to identify the origin of the n-type tendency in bulk Bi2Se3: The 
Se vacancies (Vs e i and Vs e 2) and Seei antisite dominate the donor like doping with low formation 
energy, while the predisposed Bis c i defect results in the pair of Vsei and Bi interstitial, which is 
also a donor rather than an acceptor. Moreover, for Bi2Se3(lll) surface, we find that the band 
structures modulated by the defects explicitly account for the existing experimental observations of 
n-type preference. 

PACS numbers: 71.20.-b, 71.70.Ej, 73.20.At 



Because of their excellent thermoelectric performance 
at room temperature, the layered compound Bi2X3 
(X=Se, Te) , have long been attracting great interest for 
decades in condensed- matter physics [l|, l2|- Recently, 
these materials have been drawing renewed attention by 
their novel properties as topological insulators, which 
possess insulating gaps in the bulk and gapless states on 
surfaces Specially, the defects in Bi 2 Se3 or Bi 2 Te 3 

dramatically dominate its thermoelectric and topological 
properties. For Bi2Te3, it can be simply made n or p 
type through variation in the Bi:Te ratio. For Bi2Se3, 
however, one of the major issues has been the difficulty 
in making the material p type. As revealed by angle 
resolved photoemission spectroscopy (ARPES), the as- 
grown Bi2Se_3 samples without doping are n-type semi- 
conductors [6|-|8|, mainly owing to the Se vacancies (Vg 
defects), or Se substituted defects on Bi sites (Seei an- 
tisite defects). While the acceptor- like defects such as 
Bise antisite defect, are rarely observed. Only through 
doping with extrinsic Ca on the Bi site, one can realize 
the formation of p-type material [9|. 

There have been several existing literatures dealing 
with the defects in Bi2Se3 while previous studies 

did not investigate its defect structure and formation en- 
ergy in detail. Therefore a comprehensive studying of 
the native point defect formation in this material is still 
lacking. Particularly, the influence of the native point 
defects on the topological surface states of this material 
is not much clear yet. In order to deeply understand the 
novel properties of Bi 2 Se 3 , it is instructive to investigate 
the formation of native defects, which is not only of fun- 
damental conceptual importance, but also paves the way 
for realizing promising application of Bi 2 Se3 as a novel 
topological insulator. 

In this letter, by means of first-principles calculations, 



we present a comprehensive study of native point defects 
in Bi2Se3, including Vg e vacancies, Seei and Bis c antisite 
defects. The structure and formation energy of such de- 
fects in bulk Bi2Se3, as well as their band structures, are 
calculated and discussed in detail, confirming the n-type 
preference of Bi2Se3. Besides, by analyzing the surface 
charge distribution and band structure of Bi2Se3(lll) 
surface with various defects, we present the influence of 
native point defects on the topological surface states. 

The calculations are performed within density func- 
tional theory using the Vienna ab-initio simulation pack- 
age (VASP) [13]. The PBE [ll[ generalized gradient ap- 



proximation and the projector-augmented wave potential 
[l2| are employed to describe the exchange-correlation 
energy and the electron-ion interaction, respectively. The 
spin-orbit interaction (SOI), which has been confirmed to 
play an important role in the electronic structure of the 
topological insulator, is included during the calculation. 

The crystal structure of Bi2Se3 is rhombohedral with 
the space group D^ d (R3m). Equivalently, it can be rep- 
resented in terms of a hexagonally arranged layer struc- 
ture, which is shown in Fig. [Ha). The hexagonal unit 
cell has three sets of quintuple layers, where each quin- 
tuple layer consists of five atoms, forming a stable unit 
Sel-Bi-Se2-Bi-Sel with strong intrabilayer bonds, while 
the inter-layer bonding is much weaker. Defect calcu- 
lations are performed in a 3x3x1 supercell with 135 
atoms, and integration over the Brillouin zone is done 
using the Monkhorst-Pack scheme [IH with 5x5x2 grid 
points. The cutoff energy for the plane wave expansion 
is set to 400 eV, and the structure is fully optimized un- 
til the maximum residual ionic force is below 0.02 eV/A. 
For the surface calculations, the slabs are modeled by a 
6-quintuple layer with the defects adjacent to the surface. 
The calculated structure parameters of perfect Bi 2 Se3 are 
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FIG. 1: (Color online) Atomic structure of bulk Bi2Se3 (a), 
and inter-layer structure of perfect bulk (b), Vsci (c), and 
Bisci (d). 



a = 4.198 A, c = 29.643 A, which are close to experimen- 
tal data a = 4.143 A, c = 28.636 A0. 

The formalism describing the formation energy of point 
defects in materials has been well established, which de- 
pends on the chemical potentials as well as the Fermi 
level for charged defects. The ranges of chemical poten- 
tials for Bi and Se should be evaluated as 
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where fi^i and fj,g Q are chemical potentials for Bi and Se 
in Bi2Se3, /ig; and /Ltg e are chemical potentials of pure Bi 
and Se, respectively, and AHf is the formation enthalpy 
of Bi2Se3. The formation energy of a defect X in charge 
state q is defined as [l6| 

E f [X q ] = E[X q }±fi x -E[bu\k} + q(E v + AV + e F ), (3) 

where i5[AT 9 ] is the total energy of the system containing 
the point defect X, [ix is the elemental chemical po- 
tential with a positive sign for vacancies and a negative 
sign for interstitial defects. i?[bulk] and Ey are the total 
energy and valence-band maximum (VBM) of the per- 
fect supercell, respectively. The correction term AV is 
added to align the reference potential in the defect su- 
percell with that in the perfect bulk, and Ef is the Fermi 
level referenced to Ey- Apparently the defect formation 
energy is function of the atomic chemical potential and 
Fermi level. 

Table H] summarizes the formation energy of native 
point defects in various charge state in Bi2Se3 at the 
VBM (e F =0 eV) at the Se-rich limit (/4 Se =//°J and Bi- 
rich limit (AiBi=A*Bi)' respectively. It is noticeable that for 



TABLE I: Defect formation energies (in eV) in Bi2Se3 at the 
Se-rich limit (nse=Hs e ) and Bi-rich limit (/^Bi— MbO; respec- 
tively, at VBM (e F =0 eV). 
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a specific defect the differences in formation energy be- 
tween different charge states are very tiny, which means 
that the defects in Bi2Se3 are not evidently charged. Ac- 
tually, by applying Bader analysis [TtJ , we can determine 
the ionic charges of Bi and Se atoms. The value of the 
calculated Bader charge of the individual Bi, Sel and 
Se2 atoms in perfect Bi 2 Se 3 are 0.95e, -0.55e and -0.80c, 
respectively. It is seen that the differences in electroneg- 
ativities of these atoms are not large enough to induce 
charged defect states, quite different from that in other 
typical semiconductor materials such as ZnO [HI . There- 
fore, in the following discussion we will only focus on the 
neutrally charged defect states. 

We can find from table [J that at the Se-rich limit there 
is only one stable defect, i.e., the Seei antsite that has 
negative formation energy -0.14 eV, consist with the pre- 
vious experimental observation Q. At the Bi-rich limit, 
obviously the Se vacancies with negative formation en- 
ergy are stable, and Vs G i (-0.59 eV) is more energetic fa- 
vorable than Vs e 2 (-0.23 eV). Correspondingly, stronger 
structural relaxation in the vicinity of Vg G i can be found 
compared with that of Vs C 2 ■ For Vs e i the inter-layer dis- 
tance is enlarged to be 3.02 A(see Fig. [He)) compared 
with 2.99 A(see Fig. QJb)) of perfect bulk. While surpris- 
ingly, the Bisd antisite defect at the Bi-rich limit, which 
has been considered as the p-type defect thus difficult 
to form [H, [ill , possesses of the lowest formation energy 
-1.37 eV. As illustrated in Fig. Hfd), an outwards relax- 
ation of about 1.44 A for the substituted Bi atom with 
respect to the Sel layer is observed, therefore actually 
the predisposed Big ol defect results in a pair of Vg i and 
Bi interstitial. 

In the following we present the electronic-structure re- 
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FIG. 2: Band structures of the bulk Bi2Se3: (a) for perfect 
bulk, (b) with Vsei, (c) with Seai, and (d) with Bisd- The 
Fermi level is set to zero. 

suits for bulk Bi2Se3. Figure [2] shows the electronic band 
structure with and without point defects. For the perfect 
bulk, our calculation yields a direct band gap of 0.08 eV 
at T point (see Fig. dja)), where the difference with the 
experimental value of 0.2-0.3 eV [22|, [23| may arise from 
the deficiency of pseudopotentials adopted. When de- 
fects introduced, it is clear that the gap remains robust, 
while the Fermi level of the system undergoes an evident 
shift in energy, presenting somewhat metallic behavior. 
For the bulk with defects Vs c i and SeBi, the Fermi level 
shifts up for 0.18 and 0.40 eV, respectively (see Fig. [2jb) 
and (c)). These results are consistent with previous con- 
clusions that Bi 2 Se 3 is naturally n doped [f3-|l]- While 
for the predisposed Bis i defect, the Fermi level shifts up 
as high as 0.34 eV, implying that Bi2Se3 is also n-type 
doped. Keeping in mind the much low formation energy 
and the atomic structure of this kind of defect mentioned 
above, it is much similar to that in some II- VI semicon- 
ductuors [20, [21| , the pair of Vs c i and Bi interstitial in 
Bi2Se3 forms a stable donor complex (DX) rather than 
an acceptor. 

TABLE II: Defect formation energies (in eV) in Bi 2 Se 3 (lll) 
surface for neutral charge state at the Se-rich limit (/xs e =A t s e ) 
and Bi-rich limit (/iei=/43i), respectively, at VBM (ef=0 eV). 
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0.39 
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2.57 
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Se B i 
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3.26 



Now that the basic properties of the native point de- 
fects in bulk Bi2Se3 have been established. We may ask 
how these defects affect its topological surface states, as 
the key important properties of this topological insula- 
tor. There exists a rigid difference of the chemical poten- 
tial by about 200 meV between the experimental result 




FIG. 3: (Color online) Surface charge density distributions (in 
e/A ) at the height of 2.0 A above the Bi2Se3(lll) surface: 
(a) with Vsei, (b) with Bisci, and (c) with SeBi. 

for the naturally n-type sample and theoretical calcula- 
tion for the stoichiometric Bi2Se3(lll) surface 0- The 
physic behind this difference is desirable to be investi- 
gated. Urazhdin et al. [3] have observed regular clover- 
shaped features on Bi-doped Bi2Se3(lll) surface, and 
attributed them to the antisite Bis c i at the bottom of 
the first quintuple layer. Nevertheless, we believe that 
there must be other type of defect on the surface. Based 
on the results about the defects in bulk Bi 2 Se3, we calcu- 
late the atomic and electronic structures of the surfaces 
with various defects, including Vs c i, See;, Bis c i, as well 
as the perfect Bi 2 Se3. In table [TT] we list the formation 
energy of the selected defects on Bi 2 Se3(lll) surface for 
neutral charge state. It can be seen that See; at Se-rich 
limit (0.35 eV) and V Sc i at Bi-rich limit (0.39 eV) turn 
out to be unstable compared with that in bulk Bi2Se3. 
While the defect Bisoi possesses the lowest formation en- 
ergy (-0.34 eV) at Bi-rich limit. Moreover, it is noticeable 
that for Bisd, unlike in bulk which forms a DX-like pair, 
the Bi atom takes the position of Sel atom, forming the 
surface antisite defect Big c i- 

Figure [3] presents the surface charge density distribu- 
tions of various defects near the surfaces. It is clear that 
there exists an electron density depletion at the posi- 
tion of Vsci (see Fig. [3Ja)). While for the Bis c i (see 
Fig. EJb)), the substituted Bi atom on Sel position can 
be easily identified by the electron density accumulation. 
Because of the low formation energies on (111) surface, 
and stabilities in bulk Bi2Se3, these two kinds of de- 
fects should be easily observed by STM experiment in 
Bi-doped sample. As for the Se-doped Bi2Se3(lll) sur- 
face, the Seei defect appears as a regular clover-shaped 
feature (see Fig. EIc)), with a clear depletion of electron 
density compared with the perfect surface, which can also 
be identified by experimental measurements. 

In Fig. |H we plot the calculated band structures of the 
Bi 2 Se3(lll) surface with various defects. It can be seen 
that for the perfect surface without any defect (see Fig. 
S](a)), the topological surface states form a single Dirac 
cone at the T point, implying non-trivial nature of this 
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FIG. 4: (Color online) Surface band structure of the 
Bi2Se3(lll) surface: (a) clean surface, (b) surface with Vsei, 
(c) surface with SeBi, and (d) surface with Bisei. The Fermi 
level is set to zero. 

material. Moreover, the Fermi level is -0.07 eV below the 
Dirac cone, hence keeping far away from the bulk con- 
duction band, which is much different from the recent 
ARPES experimental observation that the Fermi level 
resides in the conduction band @, exhibiting naturally 
n-type feature and presenting the existence of intrinsic 
defects. To elucidate the experimentally observed shifts 
of the Fermi level, let us investigate the band structures 
of the systems with defects. For the Vsd (see Fig. 0] 
(b)) and Seei (see Fig. H|(c)), the Fermi level of the sys- 
tems result in rigid shifts up about 0.24 and 0.10 eV com- 
pared with the perfect surface, respectively, reaching into 
the bulk conduction band, hence presenting as donors. 
Nevertheless, for Big el , the Fermi level almost overlaps 
with the Dirac point, therefore it does not present evi- 
dent electron-donor or -acceptor feature. The difference 
with that in bulk Bi2Se3 may arise from that here the 
substituted Bi atom just takes the position of Se atom, 
rather than forming the Bi-Vsei pair. 

In conclusion, we have systematically studied the prop- 
erties of native point defects in topological insulator 
Bi 2 Se3. For bulk Bi 2 Se 3 , Vs e i and Bis c i defects at the Bi- 
rich limit were identified to be energetic stable. And the 
existence of defects strongly modulate the band struc- 
tures of Bi 2 Se 3 . Moreover, we found that such defects 
have dramatically influence on the surface states of this 
material, with responsibility for the ra-type preference of 



Bi 2 Se3(lll) surface. 
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